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Abstract 

A continual model of non-singular screw dislocation lying along a straight infinitely 
long circular cylinder is investigated in the framework of translational gauge ap- 
proach with the Hilbert-Einstein gauge Lagrangian. The stress-strain constitutive 
law implies the elastic energy of isotropic continuum which includes the terms of 
second and third orders in the strain components. The Einstein-type gauge equation 
with the elastic stress tensor as a driving source is investigated perturbatively, and 
second order contribution to the stress potential of the modified screw dislocation 
is obtained. A stress-free boundary condition is imposed at the cylinder's external 
surface. A cut-off of the classical approach which excludes from consideration a 
tubular vicinity of the defect's axis is avoided, and the total stress obtained for 
the screw dislocation is valid in the whole body. An expression for the radius of 
the dislocation's core in terms of the second and third order elastic constants is 
obtained. 
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1 Introduction 



The translational gauge approach based on the Hilbert-Einstein gauge Lagrangian has 
been proposed in P for description of static dislocations in continual solids. The group of 
translations of three-dimensional space T{3) ~ IR^ is accepted in as the gauge group. 
The model pP leads, in linear approximation, to so-called modified defects instead of the 
ordinary dislocation solutions of theoretical elasticity. The modified defects demonstrate 
a non-singular behaviour, i.e., are characterized by absence of the axial singularities in- 
herent to the classical screw and edge Volterra dislocation solutions. The present paper 
is devoted to further development of the approach jlj. More specifically, it is to continue 
the investigation of the modified screw dislocation obtained in and to propose a way 
of derivation of second order corrections to its stress field. 

The point is that the Einstein- type gauge equation arising in P to govern the r(3)- 
gauge fields admits, in linear approximation, two short-ranged solutions (so-called, mod- 
ified or gauge stress potentials) which coincide asymptotically with the stress potentials 
(i.e., with the Prandtl and the Airy stress functions) of the ordinary screw and edge 
dislocations. Accordingly to the picture proposed in Ij, the stress fields calculated by 
means of the modified stress potentials just imply additional 'gauge' contributions to the 
corresponding stress fields of appropriate classical dislocations considered as background 
"configurations" (i.e., as pre-imposed sources of internal stresses). In other words, su- 
perposition of two stress fields, one is due to a chosen classical Volterra dislocation and 
another is due to the corresponding short-ranged gauge stress potential (which is local- 
ized within a vicinity of the background defect's axis), should be considered as the total 
solution of the gauge model in question. 

Therefore, two total solutions obtained in in the super-imposed form are charac- 
terized by core region, where singularities of the classical edge and screw dislocations are 
smoothed out. In other words, the gauge approach which is based on the Hilbert-Einstein 
gauge Lagrangian allows to avoid the artificial singularities of the classical elasticity. Thus, 
the gauge approach "generates" a length scale in a continuous description [T], [TT], [22] • 
The length scale characterizes the size of the domain where the classical law 1/p of the 
dislocation stresses ceases to be valid and where the axial singularity is "avoided" . Out- 
side such domain the components of the background stresses become dominating. Thus, 
in the framework of p[] it is possible to study the modified defects which allow to repro- 
duce the stresses of the classical dislocations ((T,^^ of the screw dislocation, and Upp, 
of the edge dislocation; a^^j), Ozz of the modified edge defect [T] behave unconventionally) 
sufficiently far from their axes while the stress components tend to zero within the core 
regions. 

Let us turn to the screw dislocation. Approaches j2], 0, |S] are known as attempts 
to go beyond the linear elasticity in description of the edge and screw dislocations. Non- 
hnear approach (second-order elasticity, in fact) can be used to find corrections to the 
rule However, it is still impossible to approach to the axis of a line defect sufficiently 
close since the conventional theoretical elasticity fails. For instance, the fields of second 
order stresses have been found in j2] , (by means of the stress function method) and in 
^ (in the displacement function approach) which are valid within a hollow cylinder with 
the outer radius, say, pe and the inner one pc- Free parameters of the models are fixed by 
requirements of stress-free boundaries at p = pe and p = pc- (Jpp\p=p^ = 0, app\p=p^ = 0, 
where app is the radial stress component (however, the boundary conditions are written in 
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|3] and in |1] with respect of the final and initial states, accordingly). Besides, vanishing 
of averaged over bulk's cross-section is also used for determination of one of the free 
parameters. Approaches mentioned do not consider the region < p < pe- 

A discussion of relevance of second order effects in theoretical elasticity for physics of 
imperfections in crystals, namely for modelling dislocations, can be found in j6j. Thus, it 
is clearly interesting to investigate the model proposed in pP in second order also. This 
is just the problem to be studied in the present paper. However, its purely mathematical 
aspect is of primary interest here. 

As to the gauge approaches to defects in continual solids, an attempt [7] is known to 
follow [8 in obtaining second order contributions to the stress field of the screw dislocation. 
To this purpose, the quadratic translational gauge Lagrangian |HI is used in |7j. However, 
as it is explained in pp, the quadratic r(3)-gauge Lagrangian advanced in [H], P is 
inappropriate since it forbids a modified stress potential which correctly reproduces the 
stress field of the edge dislocation. From the point of view of the Refs. ^U], ^I] ^, the 
Lagrangian used in [7], being considered as a form quadratic in the torsion components, 
is incomplete. Besides, the elastic energy is also taken in [7j in a restricted (in comparison 
with that of the classical, i.e., non-gauge approaches [2], 0, form. Since the gauge 
Lagrangian in j7j is inappropriate to capture the edge dislocation, it is also insufficient 
to consider second order corrections to the screw dislocation: the Kroner ansatz for the 
second order stresses of the screw dislocation is just of the same form as that used for an 
edge dislocation. Thus, the experience of jTj looks unsatisfactory. 

The present paper is to demonstrate that second order consideration can be carried 
out along the line of the classical investigation [Hj for the gauge model proposed in also. 
Namely, we shall consider the second order solution found in j3] for the straight screw 
dislocation lying along cylindric body as a background source of internal stresses. In 
this case, solution of the Einstein-type gauge equation gives a short-ranged "correction" 
to the classical background. The short-ranged gauge solution depends on several free 
parameters. We shall adjust these parameters in a way which differs from that in jS], jl]: 
for instance, the vanishing boundary condition will be imposed only for the outer surface 
of cylindric body containing the dislocation. Instead, we shall require vanishing of certain 
coefficients in the short distance expansion of the second order stress potential. 

Specifically, it will be demonstrated that one of the 'matching conditions' for the free 
parameters results in an expression which relates the radius of the domain of localization 
of the defect's density profile to some second and third order elastic moduli. The main 
statement of the paper thus reads: the second order solution obtained demonstrates that 
singularities in Cpp, a^tp do not appear, and the stress components tend either to zero or to 
constant values (accordingly to the choice of the free parameters) at p — > 0. However, azz 
is still weakly (i.e., logarithmically) divergent, though it is integrable over the cylinder's 
cross-section. Thus the value for cr^^ averaged over the cross-section surface is finite. 
The week divergency of cr^^ is due to a simplifying assumption about the defect's density 
profile. 

gauge approach close to ours is proposed in which is based on the translational gauge 

Lagrangian Ct written as a combination of terms quadratic in the torsion components (i.e., in the 
dislocation density's components). For a special choice of the parameters, Ct is equivalent to the Hilbert- 
Einstein Lagrangian 1 . After |12 | it is known that extension of the Hilbert-Einstein Lagrangian by terms 
quadratic in torsion (and curvature) leads to quadratic in torsion Lagrangians of more general form (as 
well as to the most general eight-parameter three dimensional Lagrangian |12|'). 
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In the present paper we are not to discuss in specific details such a comphcated field as 
description of the core structures of the crystal dislocations. Instead, only a list (inevitably 
incomplete) of further references is proposed: for instance, one should be referred to , 
|14j . [TH] . jini for the first attempts to incorporate discreteness for consideration of the 
core structures. Further references, say, for (non-linear) elasticity, for crystallography, for 
discrete (atomic) and mixed approaches, etc., can be found also in [T7j, [TH], [1^], PUI . 
For theory and experiment concerning the dislocation core structures and for effects of 
infiuence of the dislocation core structures on various physical properties of solids one 
should be referred to j2I], [22] (besides, certain refs. omitted below should be found in 

E) 

The paper is written in six sections. Section 1 is introductory (see also jl] for motives 
of our approach and for appropriate refs.). Section 2 is to outline the Einstein-type gauge 
equation. Section 3 is devoted to further specifications of the gauge equation, and a 
perturbative scheme is set up. Solution to the gauge equation which describes the modified 
stress potentials of second order is obtained in Section 4. The corresponding components 
of the stress field and their asymptotics are investigated in Section 5. Discussion in 
Section 6 concludes the paper. Details of the calculation are provided in Appendices A 
and B. Bold-faced letters are used to denote tensors of second rank (i.e., loosely speaking, 
matrices) . 

2 The Hilbert— Einstein gauge equation 

The aim of the present paper is to deduce the stress field of second order of the modified 
screw dislocation obtained, in linear approximation, in PJ. It should be reminded that 
before than in , the modified screw dislocation was already obtained in j2S| for a trans- 
lational gauge model based on the gauge Lagrangian quadratic in the torsion components 
(a 'restricted' choice of the gauge invariant quadratic form). Besides, the same modified 
screw dislocation was reproduced also in [TT], [21] for the gauge Lagrangian taken as a 
more general (in comparison with [221) quadratic form. For an appropriate choice of the 
parameters, the Lagrangian ^T], [2^ is equivalent to the Hilbert-Einstein Lagrangian 
proposed in [T|. It should be pointed out that the same, i.e., like in [2S1 and [TT], non- 
singular screw dislocation first, seemingly, appeared in [2S] in the framework of non-local 
elasticity approach. 

The main idea behind all these gauge attempts, [23], PQ? [U]? can be summarized as 
follows. Conventionally, the ordinary dislocations are characterized by the stress tensors 
cr which are singular on the defect's axes. In the gauge approaches mentioned, additional 
gauge contributions to the stress fields appear so that within compact regions (the core 
regions) the classical singularities are smoothed out. At sufficiently large distances, the 
stress components of the modified defects demonstrate the behaviour inherent to the 
classical dislocations. Within the cores transition between two asymptotics occures. 

We are going to consider the model proposed in [T], and second order elasticity ap- 
proach is accepted below. In the present section, the Einstein-type gauge equation [T] is 
outlined. Some differential-geometric notations are reminded, but for more details about 
them one should refer to [20] 5 [20]; [2Z]- It is important that now we are using the Eulerian 
picture instead of the Lagrangian one accepted in fl. The Lagrangian and the Eulerian 
pictures are indistinguishable in the linear approximation. 
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Our picture is based on the Eulerian strain tensor [20] related to deformed (final) state 
of a dislocated body. Let us denote the squared length element between two neighboring 
points before deformation as rfS"^, and the squared length in a final state will be denoted as 
ds^. We consider the difference between ds^ and dS'^, and thus we introduce the Eulerian 
strain tensor Cah as follows. 

Let us introduce the triples {x*} and {x"} as the coordinates bases (Cartesian or curvi- 
linear) to be used for description of initial and final states, respectively. The corresponding 
squared length elements can be written as gijdx'^dx^ or riabdx"'dx'' (with rjab = gijS^S^^) 
with respect to {x*} or {x'*}, accordingly. Let us define the frame components by means 
of the relation di = e"9a (here and below partial derivatives d/dx'' are denoted as 9j), 
and CO- frame components £^ - by means of the one-form dx"^ = S^dx"'. The components 

with their duals e'^j are orthogonal in the following sense: 

(throughout the paper repeated indices imply summation). Further, let us consider a map 
from an initial state to the deformed state {^''} as follows: 

Consideration of the difference 

rfs2 _ _ ^^^dCd^'' - Qijdx'dx^ = 2eabdCd^\ (2.1) 

where 

2eab = Vab - gab, Qab = gijB^B^^ , (2.2) 

allows to define the Eulerian strain tensor Cab- Here are the coefficients in 1-forms 
dx^ = B^d^"". The metric tensor gab (2.2) is called the Cauchy deformation tensor. 
In the absence of defects, B^ is expressed as follows 



B:^ — =8:-8:Vau\ (2.3.1) 

{ri) 

where = x"^ + u'^£a "with respect to the base {x*}. The covariant derivative Va in (2.3.1) 
is defined by the requirement that the components are covariantly constant, i.e., 

Va = daS; -[ah] = °' (2.3.2) 



and thus the metric rjab = ^a^bi is covariantly constant. With the help of (2.3.2), we can 
express the Christoffel symbol of second kind 1*5 1 through the metric rjab as follows: 

{ a6 } ^ ^V'''' {dar]be + dbr]ae - deTlab) ■ (2.3.3) 

In the presence of translational defects we put new Ba in another, in comparison with 
(2.3), form: 

= 1^ - < = ^; - v. + A . (2.4) 
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Here (p^ are the translational gauge potentials which are transformed under a non- 
homogeneous T(3)-gauge transformation x* — ^ x* + r]^{x) as follows: 



2.5) 



UX UX" / ■ U'lj \ 

dx^ dif 



dx^ 

The replacements (2.5) ensure the gauge invariance of the components B^. I.e., in the 
presence of defects, behave like dx^ /d^"" (2.3.1) at homogeneous [drf jdx^ = 0) trans- 
formations. 

Eventually, when (2.4) takes place, the strain Cab (2.2) can be written as follows: 

iv) iv) (iri) \ fiv) \ , ^ 

2eab = Va + ^ah + Vft Ma + ^ha " I Va ^^b + V^afe I I Vfe Ua + V'ba I , (2.6) 

where the components i^ah of the gauge potential 95 are defined by means of the represen- 
tation ip^ = £^ (f^ provided r]ab is used to rise and lower the indices. When the gauge 
potential ip is zero, Eq. (2.6) is reduced to conventinal expression for the strain e in the 
Eulerian picture pUj . 

Further, in order to consider the gauge equation, we need the Riemann-Christoffel 

a 

•ahc 



curvature tensor R - • 



^-^'-At] -Mil A'LWl] -\te\\l] . f^-^) 



where 



a 



\tc \ 

^9 (2.8) 

{6c, e}g = - {dbgce + d^gbe - degbc) ■ 



2 

In (2.7) and (2.8), the Christoffel symbols of first and second kind, {bc,e}g and 



d 
be 



accordingly, are calculated with respect to the metric gab (2.2) (the subscript 'g'). The 

(9) 

metric gab is covariantly constant with respect to the metric connection V which is ex- 

ia) 

pressed through the Christoffel symbols (2.8), i.e., equation Va gbc = is fulfilled. 

Since the metric rjab is also covariantly constant in terms of the corresponding metric 

iv) 

connection V, 

(v) , , 

Var/fec = 0, (2.9) 

we obtain from (2.9) with the help of (2.2): 

dagbc= { 9ec^{ " 1* gbe - a ebc , (2-10) 



f « 1 




f ^ 1 


I J 




[ ac j 
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where the Christoffel symbols of second kind are defined by means of (2.3.3). In its turn, 
Eq.(2.10) (plus two other equations due to cyclic permutations of the indices) gives us: 




(2.11) 

Av) iv) M \ 

"^eah = 5'''% Va ebe+ Vfe Cae- Ve Cafe j . 

We substitute (2.11) into (2.7) and obtain another representation for the curvature 
tensor: 

^ahc = Rab/ - 2 Va Cfe/ - 2eae^eb/ - (a < — >b)], (2.12) 



where Rabc'^ is the Riemann curvature calculated analogously to (2.7) but for the metric 
r]ab- In our Eulerian approach, we assume that the geometry of the deformed state is fiat, 
iv) , 

and so we put Rabc equal to zero. Now we are ready to write the Hilbert-Einstein gauge 
equation. Let us define the Einstein tensor G"^-^ as follows: 



Qef ^ ^Seabgfcd-^^^^^^ (2.13) 



or, with the help of (2.12), 



iv) iv) 

Qef ^ _geabgfcd _ 2S'^'S^'''eade' e,/ . (2.14) 

In (2.13) and (2.14), E'^^'^ is the totally antisymmetric Levi-Civita tensor [20] defined by 
means of the metric rjab- Therefore, the gauge equation proposed in ^ (Section 6) takes 
the form: 



('?) iv) 

_geabgfed ^^^^ ^ {2s)-\a''^ - (^bg)*^^) + 2S'''''£f''''eade' Cb/ . (2.16) 

Variational derivation of (2.15) can be discussed along the line of '1] where the gauge 
approach was developed in the Lagrangian coordinates. Right-hand side of (2.15) is given 
by the difference cr — cTbg, where crt,g implies the stress tensor of a background defect. The 
difference cr — cTbg, i.e., just the deviation of the total stress from cTbg, plays the role of the 
source of geometric configurations described by the Einstein tensor G. The parameter 
s (a coupling constant, accordingly to gauge terminology) characterizes an energy scale 
intrinsic to the gauge field (p, and it appears as a factor at the Hilbert-Einstein Lagrangian 
density p. 

In the present paper, cTbg is assumed to be given by the stress field of a single straight 
screw dislocation lying along an infinitely long cylindric body. Practically, the solution 
provided by [3, will be adopted, which is valid within a hollow cylinder restricted by two 
surfaces: internal (p = pc) and external (p = pe). Further comments about (2.15), (2.16), 
and about their specification for the present problem can be found below. Besides, both 
<T and CTbg are assumed to respect the equilibrium equations: 

iv) , iv) , 

Vaa'^' = 0, Va (o-bg)'^'' = 0. (2.17) 
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More detailed consideration of the gauge geometry behind the model to be studied in 
the present paper should be done elsewhere. However, Refs. ^28, should be mentioned in 
addition to those listed in jT] since they contain reviewing notes and useful refs. concerning 
translational gauge geometry. It is also interesting to note Ref. [221 where a topological 
picture is build up which includes dislocations (torsion) and exta- matter (non-metricity) . 



3 Specification of the gauge equation 
3.1 The stress function method 

We shall investigate Eq. (2.16) using the method of stress functions proposed in [2] for 
solving the internal stress problems in incompatible elasticity. Specifically, we shall follow 
Ref. 0, where the approach ^ was further developed in the successive approximation 
form to determine the second order stress fields of the screw and edge straight dislocations 
lying along cylindrical tubes of circular cross sections. An exposition of [H] can be found 
in Ref. [20] devoted to a review of dislocation problems in non-linear elasticity. Certain 
details (concerning, for instance, the tensor formalism in curvilinear coordinates) omitted 
in what follows can be restored with the help of [3] and j^Oj. In what follows, bold-faced 
letters denote tensors, and all the indices can easily be restored. 

Accordingly to Eq. (2.6), let us represent the strain and the stress tensors as the 
perturbative expressions: 

(1) (2) (1) (2) ,^ , , 

e = e + e , cr = cr + a , (3.1) 

(2) (2) (1) (1) 

where e , cr are assumed to be of second order smallness in comparison to e , a . Ex- 
pressions (3.1) can be understood as the first two terms of formal perturbative series in 
powers of a small parameter. Substituting (3.1) into (2.16) and (2.17), we obtain the 
following two sets of the governing equations: 



(first order) 



(1) 

v„ V '^^ = , 



(second order) 



(inc 



Va = , 



(3.2) 



/ (2)\ ab , , 1 / (2)\ ab W , 

(inc e) = (2s) -1(5 a) + Q '''' 
The following notations are used in (3.2) and (3.3) {i = 1,2): 



(3.3) 



(inc ^e^)"' = -S'^'^'^S'^'V^Vf e^, , (3.4.1) 
{6^y'^'S^''-(^,,y\ (3.4.2) 
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K J uf (1) (1) ; , . 

where 'Inc ' denotes, so-called, incompatibility operator (acting on a tensor argument), S"''"^ 

iv) 

is the Levi-Civita tensor, and Va implies the covariant derivative Va- Indices in (3.2)- 
(3.4) are raised and lowered by means of the metric r]ab- A direct comparison, for instance, 
with Eqs. (622) and (623) provided by [20] demonstrates that the second equations in 

each pair (3.2) and (3.3) imply transparent modification of the corresponding equations of 

(i) 

the conventional approach: the only difference is given by the terms (2s) 6 cr^ which are 
responsible for the short-ranged behaviour of the resulting stress functions. Besides, the 
contributions due to the contortion (and thus due to the dislocation density) are absent 
(see explanations in [TJ). 

The elastic energy potential is chosen in the Eulerian representation as follows 0: 

W{e) = jll{e) + kh{e) + Vll{e) + m'h{e) h{e) + n'h{e) , (3.5) 

where j = fi + A/2, k = — 2/i (A and /i are the Lame constants), while l',m',n' are the 

elastic moduli of third order. For the given choice of the potential W{e), the constitutive 

(i) (i) 

law which relates e to cr takes the following form [21 : 

(1) „ ^ „ (1) 
e = Cili[cr)r] + C4 cr , 

(3.6) 

(2) ^ J .(2)^ , ^ (2) , « 

e = Cili{cr)r] + C4 cr + ^ , 



where 

^= [C^l'tC^') + C3/2(iT')jT7 + C^hC^) a +C,hXcr)[^^) \ (3.7) 
and the numerical coefficients are 

C2 = -u)-l) + 3L + M, 

C,^^^{l-u') + M, (3.8) 
C5 = -^(3z.-2)(l + z/)-M, 

In (3.8) we are using the Poisson ratio u and the elastic modulus E: 

A p /i(3A + 2/i) 



2(A + /i)' A + /i ' 

besides, the relationship between the elastic parameters of third order L, M, N in (3.8) 
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and m', n' in (3.5) is given as follows: 
n' + S/i^A^ = 6/i , 

3m' + n' + 12/i2/C(3M + iV) = , (3-9) 

1 + 

27/' + 9m' + n' + 27/C^(27L + 9M + A^) = 36/C , 

where 

2 

3(1^ = ^+3^- 

Besides, 77 implies the metric rjab in (3.5)-(3.7), and the functions Im {rn = 1,2,3) of a 
tensor argument, say, t are defined as follows: 

/i(t) = tr(t), 

m ^ Uim - hit')) , (3-10) 



2 

/3(t) = Det(t). 

The Cayley-Hamilton theorem must be used to express the inverse (ir) in (3.7). More 
details about derivation of Eqs. (3.6)-(3.9) can be found in jS], pUj . 

Now we are in position to use the first and the second equations (3.6) in (3.2) and 
(3.3), respectively. We obtain: 



(first order) 



(1) ab 



A + (1 - a)(V,V, - VabA)h{l) = ■ 



(3.11) 



(second order) 



Va a = ^ 



A ^a' + (1 - a){VaVb - VabA)h{2) (3.12) 



''(5 o-y + 2/x Q - 2^(lnc ^ 



(1) (i) 
where is given by (3.7), h{i) = Ii{(t), A = VaV^ and Q ^'"'> is expressed by means of 

(3.4.3) and (3.6). Besides, we make use of the parameters n"^ = fi/s and a [H], 0: 

A 1 ^ 2(A + /i) 1 

— 1 — a 



3A + 2/i 1 + ' 3A + 2/i 1 + z/ 

The curly brackets around the indices imply symmetrization. 
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We use the stress function ansatz to fulfil the equilibrium equations for stresses given 
by (3.11) and (3.12) as follows: 

<T= Inc X . (3.13) 
Substituting (3.13) into the second equations in (3.11), (3.12), we obtain a couple of 
equations to determine the stress potentials X: 



AA Xab +a(VaV6 - r]abA)Ah(x) + ((1 - a)VaVfe + ar/a6A)v"V'^ Xcd- 
- A(VaVe X \ + VfeVe X \) = K^{S a-)^^ + 2/i § , ^ = 1, 2 , 



(3.14) 



where 

(1) (2) (1) / (1) 

5" [ab) = , S (ab) = Q (ab) - (InC ^ 

(i) 

and 5 cr is written by means of (3.4.3) and (3.13). 

Using a linear transformation to another stress potential X ', 

x = x' + Yzr^ Vabh{x '), 

we can reduce (3.14) at z = 2 to a more simple form: 

AAx'-n^Sa' = 2fiS, (3.15) 

(2) (2) 

where prime at cr implies that the tensor is expressed through X ■ Provided a tensor- 
valued Green's function of the corresponding operator acting in L.H.S. of (3.15) is known, 
solution of (3.15) can be obtained in a standard way. However, in what follows we shall 
be concerned with (3.14) itself. Therefore, the main task below is to adjust (3.14) to the 
special case in question, i.e., to the case of the screw dislocation along a cylindric body 
of circular cross- sect ion. 



3.2 The gauge equations in the first and second orders. The 
choice of the model 

Owing to the fact that the equilibrium equations given by (3.11) and (3.12) (the first ones 
in pairs) are fulfilled by the Kroner ansatz (3.13), the problem is to determine the stress 

function components Xab from the gauge equations given by (3.11) and (3.12) (the second 
ones in pairs). 

Now we replace all the derivatives Va by the partial derivatives da = d/dx°', where x"" 
are the Cartesian coordinates in the final state [20j. We assume also dz = = 0. Let 
us introduce the following notations for those components of the stress potential which 
are non-trivial: 

(i) (i) (i) 

fi (j)= 82 Xi3 -di X23, « = 1, 2 , 

(3.16) 

/ = X33, P = -5n X22 -dl2 Xii +2dl2 X12 . 
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The other components Xab are zero. The background stress tensor cTbg (see (2.17)) is 
also given by (3.13), though the corresponding stress functions are labeled by appropriate 

1 ■ 

subscript: Xbg- 



The Grst order 



Since the background stress field is assumed to correspond to that of the screw dis- 
location in the form provided by we conclude that only the component of the stress 
(1) 

potential is nonzero in the first order. It is described by second equation in (3.11), and 
the latter acquires the following form PQ: 

/ (1) (1) .\ 

da[A ^ -k\(P - = 0, a = 1,2, 

or 

A = - 0bg), (3.17) 



(1) 

where 0bg= ("V^^r) log p. It is appropriate to re-express (3.17) as follows: 

(a - ^'){<P - 0bg) = b6^'\x). (3.18) 
Solution to (3.18) describes the modified screw dislocation, and it is given by 

(1) (1) 



bg 



-fs, fs = {b/2n)KoiKp) . (3.19.1) 



From (3.13) and (3.16) we obtain the only non-trivial component of the total stress as 
follows: 

cT^z = - dJfi (t>) = ^ p'\l - KpK,{Kp)) . (3.19.2) 



27r 

Equations (3.19) witness about existence of a core region at p< 1/k: outside this region 
the gauge correction to the classical long-ranged law 1/p is exponentially small. At p <C 
K,^^, the law 1/p (characterizing the stresses) is replaced by another non-singular one. 
More detailed information concerning the numerical behaviour of the solution (3.19.2) 
(including a treatment of k,~^ in terms of interatomic spacing) can be found in and 

I2n!. 

The second order 

In order to obtain the gauge equations in this case, it is necessary to specialize the 

(2) (1) / (IX 

source Sab in (3.14), i.e., its terms Q (at) and line ^] , which depend on the first 

^ ' V / (ab) 

order solution (3.19.1), have to be written explicitly. (In what follows we shall use Q, 
e without the superscript ^^\) 
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Let us begin with Q(^ab) (3.4.3). More explanations can be found in [3], [20]. Using 
(2.11) we obtain from (3.4.3) the following non-zero contributions: 



Q33 



Q 



22 



(^1623 - ^2613)' 



4 (i9ie23 ^2613 - 82623 dicis) + ((9ie23 - c?2ei3)^ , 



(3.20) 



while Q12, Q23, Q13 are zero. Our expressions (3.20) differ from the analogous quantities 
in [3j, 1201 (^ different numerical factor in Qn, Q22, and absence of the second term in Q33) 
since in our case the corresponding Einstein tensor in L.H.S. of the gauge equation (2.15) 
is expressed by means of the Riemann-Christoffel tensor only pp. However, in [201 
the corresponding incompatibility equation includes a contribution from the contortion 
tensor also. Now let us obtain ^ (3.7) in components: 



(l//^')^ii 



(1//U')^22 



(l//i')^12 



-C3 
-C3 



, (l)s2 f (l)x2 

(di 0) + (^2 0) 

^ (^\2 . (IK 2 



C7 5i 



(IK 2 



^7(^2 0) 



(1) (1) 



(3.21) 



(l/^2)v[/ 



33 



(1/^2)^21 = -Ct^i d2 
-C3 



(9i 0) + (^2 0) 



Now we are ready to consider the gauge equations of second order given by (3.14). We 
obtain: 





.(1 


- a)p + 


aA/ - 


n'if- 


/bg 




a 


- a)p + 


aAf - 




/bg 


dl2 


a 


- a)p + 


aAf - 


^M- 


/bg 


1 - a) AAf - 


f aAp 


- K^ip 


-Pbg) 





2/i Qii + 



33 



2/i Q22 + 5?i^33 



(3.22.1) 



-2/i VI/ 



33 



2/i Q33 - 5?2(^12 + ^21) + 5?1^22 + 9|2^11 



(2) 



(3.22.2) 



Besides, there exists a couple of equations to determine 0. But since \E'i3, \E'3i, ^23, ^32 

(2) (2) 

are zero, and 0bg is also zero we put consistently 0=0. Thus, the only equations to 
be considered are given by (3.22), where the corresponding components of Q and ^ are 
given by (3.20) and (3.21), respectively. 

Equations (3.22) look very similiar to those of the classical approach [3], [201 excepting 
of the fact that Qn, Q22 are zero in the classical consideration. The point here is as follows: 
equations for the classical stress potentials of second order, pbg and /bg, are considered 
for a two-dimensional domain Pc < P < Pe (Pe and pc are the external and internal radii, 
accordingly). However, Qn and Q22 are equal to 4(9ie23 — 02613)^ in [3], and therefore 
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(1) 

they are proportional to A ^^^g, while the latter represents the dislocation density profile 
(i.e., the dislocation density component ^^12)- Outside the disc given hj < p < pc, 

(1) 

A 0bg is zero because of its proportionality to the Dirac 5-function. It is just the case, 
why Qii and Q22 drop out of the classical version of Eqs. (3.22). 
Accordingly to jS], a constant Q' can be introduced as follows: 

Qn = dl,Q', Q22 = dl,Q', Qx2 = -d\^Q\ (3.23) 

where Q' is to be adjusted in the end of calculation of the stress components ^. Therefore, 
equations which define pbg and /bg take in our notations the form (compare, for instance, 
with (3.22)): 

AA/bg = 2/i^ A(vi/^f + gO + [g^f + dl,^\l + dl^^^f, - 2d\2^% , 

(1 - a) pbg + a A/bg = -2/i (v|>^f + Q') . 

(3.24) 

The first and the second equations in (3.24) correspond, respectively, to Eqs. (657) and 

(2) 

(655) in [2ni, provided (733 and F(2) therein are identified as pbg and — /bg- Besides, Eq. 
(3.23) defines Q' in opposite way (the sign is different) in comparison with Eq. (658) in 

In the present paper, Qn, Q22, Q33 are given by (3.20), and we find: 

Qii = Q22 = (^-A 0j . (3.25) 

Is it possible to find a "potential" ^ = ^(p) which is analogous to Q' of the classical 
approach? For such g = g{p) the following equations should be respected: 

Qii = 9^29 = ^{dlp + ^dp) g - ^^^^ (dip - -^dp)g , 

Q22 = dl,~g = \{dlp + -^dp) g + ^ {d% - i^,)^, (3.26) 

Qu = -dl29 = -'^{d%--d,yg=0. 

z p 

Equations (3.26) lead us to the following pair of equations: 

K-ld,)g = 0, 

{d% + hp)g = 2Qn = 2Q22. 
However, Eqs. (3.27) are not consistent for the given Qn, Q22 (3.25). 



(3.27) 



(2) 

^More precisely, the numerical value for Q' arises from a condition that a mean value of <j 33 is zero 

m- 
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(1) 

In the classical approach, A ^^^g corresponds to the defect's density profile, and the 
latter is given by the Dirac 5-function. Therefore, Qu = Q22 = at p > pc, and so it 

(1) 

is possible to choose Q' as a constant. In our approach, A also represents the density 
profile of the modified defect, and it is given as follows [T], [TT], ^B] : 



(1) 



b ^2 
— k^Ko{kp) 



(3.28) 



In the limit k —>■ 00, the function Kq{kp) demonstrates a behaviour of a 5-like function 
on a plane "centered" at p = 0. Therefore, the following simplification can seemingly 
be made to keep the situation tractable in the framework of the plane problem: let 
us approximate n"^ Kq{kp) by a piecewisely constant function which takes two different 
constant values either within the disc < p < pc, or outside it. Besides, we shall assume 
that differentiations of this "hat" -function at p = p^iO are negligible. Then, it turns out 
that equations (3.27) become consistent. Clearly, such approximation is rather rough at 
p <C K~^, i.e., in a very close vicinity of the classical defect's axis. However, as we shall see 
below, this simplification leads to a reasonable picture for a non-singular modified screw 
dislocation in the second order also. 

The replacement proposed for the density profile is given as: 



^k^Ko{kp) 
Ztt 



T^Pc 



(3.29) 



where /i[o,pd(p) is equal to unity at < p < pc, and to zero otherwise. With the density 
profile given by (3.29), we find g which respects (3.27): 



c . 



2'npl 



(pV2) + C" 



Pc< P, 

< P < Pc , 



(3.30) 



where the constants C and C" will be adjusted later. 
Let us write Qsa (3.20) explicitly: 



Q 



33 



(1) 



(1) 



d'22 



(1) 



(1)\2 



(3.31) 



Using (3.21) and (3.31), we calculate: 

Q33 - 2df,^,2 + df,^22 + dl,^u 



A^33 + (1 - 2p'C7 



(9; 



(1) 



9 .9 

2 



(ly 



22 



+ 



1 2 

-A : 

2 



(3.32) 



Therefore, taking into account (3.26) and (3.32), we obtain from (3.22) the following 
couple of equations to determine / and p: 



'l-a)p + aAf - - /bg) = -2p (^33 + 



(3.33) 



;i - a) AAf + aAp - n^p - pbg) 



p 
2 



A 



(l)x2 



+ 
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+ 2^A*33 + 2^(1-2^2(^7) 



/ . (i)x2 „ (1) „ (ly 



We obtain p from (3.33) and substitute it into (3.34): 



a 



P = 



A/ + 



K 



if - /bg) 



2/x 

1 - a ' 1 - a ^ 1-a 
(A-«:^) (A + ^)(/-/bg) =7e, 



(*33 + g) , 



where 



l-2a 
2// 



7^ = (A-«:2)(^33+^_^bg_g,)^ 

+ (l-a)(l-2/.2Cv)($-$bg)-^(A0y, 
$ = 0) - 9n ^22 ■ 



(3.34) 

(3.35) 
(3.36) 



(3.37) 



(3.38) 



Here /bg respects the first equation in (3.24), and $bg) ^af are given by (3.38), (3.21) 

(1) (1) " (1) 

(provided (f) is replaced by 0bg)> correspondingly. Besides, we formally keep (A 0) as 

an exact expression. 



3.3 Final remarks about the gauge equations 

To conclude Section 3, we should pay attention also to the structure of the function TZ 
in the righ-hand side of (3.36). Explicit expression for TZ is given by Eq. (3.37). Our 
problem has an axial symmetry, and therefore the use of the cylindric coordinates p, (p, z 
instead of the coordinates in final state x°-,a — 1,2,3 (the coordinates p and ip are chosen 
in (x^, a;2)-plane and z = x^) is more appropriate. Therefore, we obtain from (3.21) the 
following expressions in the cylindrical coordinates: 



(1), 



33 



-C{dp (f)) , C = Cg , 



A*33 = -2c 
where $ (3.38) is re-written as: 

1 



(1) 



(1) 



(1) 



d, d,{A ct>) + (A ct>y 



- 4c$ 



(1) , (1) 1 (1) (1) (a/iV 
-9, a, 0A0 +^^^ 

p p p^ 



Further, we use (3.39), (3.40) in (3.37), and obtain in the following form: 

2 / (1) 



n 



k 



W{p) + (-^ + k'c) ((5, 0)' - {d, 0bg)') 



- -(a (/> 9, - a 0bga, 0bg) 

- 2c 9p(A 0)0', 4> -9,(A 0bg)9, 0bg 



(3.39) 



(3.40) 



(3.41) 



16 



where 



2/i 



l-2a 
1 + a 



a)(l 



4c, 

(1) 



W{p) ^ (A 0) - 2c((A cf>) - (A j - K'{g - Q') 



(3.42) 



The contribution most interesting for us is given by 

I.e. 



Let us have a look at TZ (3.41 

the second term in it. The other terms in (3.41) contain either A or dp{ls. 
are dependent on the density profile or on its derivatives. In other words, these terms are 
more significant either within the core or near its boundary. 

Let us rewrite TZ (3.41) once again using the explicit expression for the background 



(1) 



solution of the first order 0bg= ("V^^r) log p. We obtain: 



n 



k 



(1) 



W{p) + + «'c) d,fs [dpfs - 2d, 



+ 



cA ^ + 2cd,{A cp)d,fs 



(1) 



(3.43) 



2c (a 



(1) 



(1) , (1) 



bg 



Now the structure of TZ can be characterized as follows. The first term, W{p), is deter- 
mined by the density profile, and it seems to be significant rather within the core since the 
constant Q' can be removed outside < p < pc by an appropriate choice of C in g (3.30). 
The second term in (3.43) behaves as at p <^ 1, and therefore it is just responsible 
for the fact that the stress function to be found / is expected to cancel exactly the most 
important term (oc log^p) in /bg at up <^ 1: 



AA(/ - /bj 



J^kc 
'4^p4 



AA/, 



bg 



kc 
4^/ 



while /bg itself is given by Eq. (5.12) below. The third and the fourth terms are concerned 
with the density profile and with its derivatives. It can be assumed that the last two 
terms imply an increment of the corresponding quantities in the gauge approach due to 
a replacement of the density profile of the background defect by the density profile of 
the modified defect. In what follows, we shall neglect a possible infiuence of the last two 
contributions. Thus we obtain: 



TZ ^ k 



W{p) + (4 + '^'c) d,fs {d,fs - 2d, 0^bg) 



(1) dJs / 
+ cA ^ + 2ca,(A (j))dpfs 



where we put approximately 

W{p) 



(3.44) 



(3.45) 
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4 Solution of the gauge equation 
4.1 Preparation 

Now the task is to solve Eq. (3.36) with the R.H.S. given by 7^ (3.44) with W{p) (3.45). 
We shall do it in two steps. As the first step we shall solve equation 



z + z — — z 

dz'^ dz 



G{z) = kn{z) 



(4.1) 



where the variable z implies the radial coordinate p rescaled as follows: z = up, d/dz = 
K d/dp. Equation (4.1) is a non-homogeneous Bessel equation [HU], and its new source 
(after multiplication of equation by p^), though again denoted by 7^, is written in the new 
variables: 

n{z) = w^{z)z^ + X^{c + cz^)Ki{z) {Ki{z) - -) 

f d \ ' ^^-^^ 

- X^zK,{z) fcA,0 + 2cz— (A,0) j , 

jj2 d'^ -1 d ~ 

where = — -, = -r-^ + z —, Ki(z) is the modified Bessel function EOj, and 

All'' dz^ dz 

(1) h 
implies with removed factor — . 

27r 

In order to explain the notation Wc{z) in (4.2), let us have a look at W{p) (3.45). Let 
us choose the constant C, which appears in g (3.30), equal to the constant Q' defined by 
(3.23). Then we obtain for W{p): 



W{p) 



h \2 1 

2^J l4 



;i + a) - y(p^ - pI) 



](P), 



(4.3) 



where C" (see g (3.30)) is fixed by requirement of continuity of the density profile at pc- 

62 1 



C" + 



Itt"^ p1 



Q'. 



After the replacement of Kp by z, we obtain for k '^W{p): 



(4.4) 



where Zc = upc, and h[o^z^]{z) is equal to 1 at 2; G [0, Zc] or to zero, otherwise. 

Solution of (4.1) is based on a knowledge of the asymptotical behaviour of 71 (4.2) 
at z ^ 1 and z 1. Let us obtain the corresponding expansions. First of all, let us 
take into account that A^^ is exponentially small at large z (or even zero, provided the 
replacement (3.29) is made) since it is equal to —Kq{z). Therefore, TZ is localized at large 
z. 

In order to study the case 2; <^ 1, let us represent TZ as a sum: 



7^ = wc{z)z^ + 7^l(2) + 7^2(^) + n^iz) . 
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We obtain the following expansions for TZi, 7^2, T^s at small z: 
|i ^ ( c + cz') K,{z) {k^{z) - I) 

- - - c+ ^zHog^z - ^(l-2log^)zHogz (4.5) 
+ -^(l-21og^)%2 + oizHog'z); 



^ = -czKi{z)l^^4) ^zKi{z)h[o^;,^]{z) 



c 

^2 



'^''hio,.M (l + l^'^ogz - i-(l-21og^) + oiz') 



4 



(4.6) 



-2cz'K,{z)^(A,cP) 



(4.7) 



'dz 

~ -2c - 2c;2^1og;2 + c(l - 2 log + o{z^) 

In (4.6), we took into account the replacement (3.29). For a comparison, the exact 
expansion for TZ2 looks as follows: 

|| -c\og{lz) + \z- (l - 2 + o{z^) . (4.8) 

It is seen from (4.8) that the terms oc log 2; and oc z"^ log^ z are absent in the approximate 
expression (4.6). 

Expansions (4.4)-(4.8) suggest the following series form of TZ (4.2) at 2; <C 1: 

n{z) ~ PoZ^ + Pl + P2 log 2; + P3 

(4.9) 

+ P4 Z^ log^ 2; + P5 2;^ log 2; + P6 + 0{Z^) , 

where po = P^iA implies wjyz) (4.4). The coefficients pi, p2, Ps, P4, Phi Pe are influenced 
by our assumptions, and we obtain them as follows: 

Pl = -X^c, 

2 ^ 

P4 = ^ ^ , (4.10) 



'2c 



— c 



" 16 



Practically, only TZi (4.5) and 7^2 (4-6) are taken into account in order to assign the specific 
values (4.10) to the corresponding coefficients in the expansion (4.9). The contribution 
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7^.3 (4.7) is excluded from the consideration since we use (3.29) and neglect possible 
contributions which should be important near the boundary of the core. However, a 
possible implication of TZ^ for the coefficients in (4.9) would be given by shifts in P3, p^, pe 
by numerical constants dependent on c (3.39) (for instance, — c in p^ would be replaced by 
—3c). A usage of TZ2 in the form (4.8) instead of the approximate expression (4.6) would 
remove in ps, ps, p^ (4.10) the dependence on the value of the core radius Zc. Besides, 
usage of (4.8) would lead to non-zero contributions in p2 and ^4. But since our choice of 
the density profile is given by (3.29), the coeficient p2 is simply zero. By formal reasons, 
we keep the corresponding terms in (4.9) with unspecified p2 which is to be equated to 
zero at the very end of the calculation. 



4.2 The solution 

Therefore, let us first consider the non-homogeneous equation: 

rf2 d 



^2_^ _|_ _ 

dz'^ dz 



viz) = kUiz) 



(4.11) 



General solution to (4.11) is given by the standard formula [HU] : 

k~'y{z) = Ayiiz) + By2{z) 

^ yi{z)y2{t) - y2{z)yi{t) n{t) 
yi{t)y'^{t) - y2{t)y[{t) 



dt. 



(4.12) 
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Here zq is to be adjusted, and yi{z), y2{z) are linearly independent solutions of the 
corresponding homogeneous equation. 

In our case we put yi{z) = Io{z) and y2{z) = Kq{z) [10], and (4.12) gives us solution 
to (4.11) as follows. Let us define G{z,s): 



k''Giz,s) = loiz) (^A{s)+ jKoit)nit) 
+ Ko{z) (sis) - J /o(^)7^(t) 



dt 

T 

dt' 

T 



(4.13.1) 



where 



dt 
T 



A{s) = - J Ko{t)n{t) 

s 

1 ^ 
const + / + y log^ + + ^3)- 



(4.13.2) 



B{s) ^ 

Then, solution to (4.11) appears as 

G(z) = lim G(z,s) = G(z,0). 



dt 



(4.14) 



Asymptotical behaviour of G{z) can be obtained from (4.13), (4.14). At large z, 
G{z) decays exponentially, i.e., G{z) is well localized. Using the expansions provided in 
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(4.15) 



Appendix A, it is straightforward to establish the behaviour of G{z) at small z: 
k ^ G{z) ~ go + + ?2 iog^ z + log^;2 + 54 log 2; 

+ ^5 z'^ log^ 2; + ge log^ z + q^z"^ log 2; + gg -^^ 
where 

go = -const X log - - - — , 

Pi P2 Pi , P3 , , 3pi 

9i = ^ , 92 = , ^3 = + y , 54 = - const + — , 

P2 Pi , P3-P2 , P4 

const pi 3p2 P3 P4 , P5 PO - , 

const / 7n 5pi p2 , 3p3 8^4 , Pe - Pb 

In (4.16), Pi, ■ ■ ■ ,P6 are given by (4.10), 2k is given by (A8) in Appendix A, and const is 
introduced by the definition of B{s) (4.13.2). It is seen that the term p2logz in (4.9) is 
just responsible for the highest powers of the logarithm in (4.15): g2log^2; and q^z^ log^ z 
(and similarly for the logarithms at higher powers of z^). 

As a second step, wc arc going to find the modified stress potential of second order 
/ = /bg + JF, where respects the Bcsscl equation 



Z^-r^ + Z— + Z^ 

az"^ az 



m = ^.G{^.). (4.17) 



In (4.17), the variable z is defined differently, z = J\fp, J\f^ = — , and G{z) is given 

by (4.13), (4.14). Using again (4.12), we obtain solution to (4.17) in the following form: 

^ y (4-18) 

/^(p) = M^fp) I Yo{J\ft)G{nt)tdt - Yo{Np) I Jo{Aft)G{nt)tdt, 
p p 

where ^0(2;) = {n/2)Yo{z), and Yo{z), Jo{z) are the Bessel functions which solve the 
homogeneous version of (4.17). 

For our purposes it is appropriate to put D — and C 7^ in (4.18). Now let us write 
the asymptotical results for J-'{p). At large p, the contribution given by Ijr is exponentially 
small, and thus YQ{J\fp) dominates in J-{p) at J\fp ^ 1: 

For small J\fp, we obtain (Appendix B): 
T{p) ~ ro + ri logp + r2 log^p 

(4.19) 

+ rsp^log^p + tap^Xo^ p + rsp^logp + tq p^ + rjp'^log^ p, 
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where 



Xy + {C + Xj) 



n = C + Xj 



r2 



Pi 



r3 



P2 

24 ' 



r4 = J5 - n + 3 log AT (J4 - n) 

= 4fl-^)--(l-log-)- + - 
[V J 32 ^ ^ ' 8 8 

= Je-% + 2\ogAf{J, - %) + 31og2Ar(J4 - %) - C 

r, = Jr-Y-r + \og^ {J, - %) + \og^ M [J, - %) 



(4.20) 



n = Jg - Yg = k 



K 



384 



-P2- 



In (4.20), Xy, Xj, and Jj, (z = 4, . . . , 9) are defined in Appendix B, and pi, 
given by Eqs. (4.10). 



,Pg are 



5 The stress tensor 

5.1 The components app and 

In the previous section we found J-'{p) which respects the inhomogeneous gauge equation 
(3.36) with the source term 7l{p) taken in the approximate form (3.44). The function J^{p) 
implies the difference /(p) — /bg(p), where /(p) is the modified stress potential of second 
order, and /bg(p) is the background stress potential. Practically, J^{p) is given by the set 
of integral representations (4.13), (4.14), and (4.18). However, we shall not attempt to 
elaborate a single formula which would express J-'{p) through TZ{p) more explicitly. The 
most important for us asymptotical properties of the stress field of the modified screw 
dislocation can be obtained just from Eqs. (4.13), (4.18). 

Equation (4.18) allows to express the modified stress potential / as follows: 



/b 



/ = /bg + CYgiAfp) + 

^^S^ P + diP^ + d2 logp. 



(5.1.1) 
(5.1.2) 



where the classical second order stress potential /bg (5.1.2) is written in the form suggested 
by Free parameters di and d2 in (5.1.2) are determined in [Sj from the requirement 

(2) 

that the second order stress a pp vanishes at the boundaries of a hollow cylinder p = pc 
and p = pe > Pc- In what follows, it will be seen that the cut-off at p = pc disappears in 
the gauge model proposed. 
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The second order stress tensor of the modified screw dislocations is given by the 
following relations: 

(2) Id (2) _ ^ (2) 



^PP^-Z~„J^ ^'t>4>^~T^J^ (^zz^P, (5.2. IJ 



(2) 

where / is given by (5.1) and p is given by (3.35). Other components of cr are zero. It 
is most important for us to consider the limiting behaviour of the solution (5.1), (5.2) at 
N'p ^ 1 and J\fp <C 1. Let us begin with the case 

A) A^p < 1 

(2) (2) 

Since our attention is attracted now to a pp and cr (p^,, let us write, using (4.19) and 
(5.1), their expansions as follows: 

- (ri + d2 — 

p^ p^ 

2(r4 + rs) log p - (rg + 2rQ + 2di) - 4rrp'^ log^ p , 



= - 2{r2- k^)^-^-(n + d,)\ -2rs log' p - (Sr^ + 2n) log' p 



(2) 



- 2r3log^p - (9r3 + 2r4)log^p - 2(3r3 + 3r4 + r5)logp 

- (3r5 + 2r4 + 2re + 2di) - 12r7pMog^p. 

(5.2.2) 

It was noticed in Subsection 4.1 that for our choice of the density profile, the series 
expansion TZ (4.9) at small distances is missing the term corresponding to the coefficient 
P2- In its turn, the vanishing of p2 just implies absence of the terms depending on and 
r7 in (4.19) and so in (5.2.2). Besides, the constant term tq is irrelevant for the stress 

(2) (2) 

components a pp, a (5.2.1). The contribution corresponding to r2 is compensated 
exactly by the first term in /^g (5.1.2) (and thus the contribution oc p^^logp disappears 
in (5.2.2)). 

Therefore, our attention should be paid only to the coefficients containing ri, r4, rg, re 
in the final expansions (5.2.2). First of all, it is interesting that there exists an opportunity 
to make equal to zero: r4 = 0. Indeed, using in the form (4.20) (where we put P2 = 
and use pi, pa in the form given by (4.10)), we can re-express equation r4 = as follows: 



or, after the use of c in the form (3.42), 

l + v i V2(l-^^) z, 



c=(^^^-4c) (,77^ + ^ . (5.4) 



c , 



Let us rewrite (5.4) using c (3.39) as follows: 
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where a = u). Left-hand side of (5.5) is positive (it is known that < < 1/2 

for isotropic materials), and thus 77 in its R.H.S. is also positive. Therefore, the following 
two requirements appear: 

a) C3 > 0, C7 < l/2/i2, 

h) C3 < 0, C7 > 1/2//^ 
Provided these requirements are fulfilled, Eq. (5.5) can be solved for Z: 

Z = (5.6) 

In its turn, Eq. (5.6) results in the following restrictions on the parameters C3, Cf. 



a) C7 + 2(l + z/)C3 < < C7 + ^^^til)!c3, 

(1 + .)^ 1'" ^ (^■^) 



h) C7 + ^ ^ < — < C7 + 2(l + i/)C3. 

(Notice that (1 + i/)/i/ > 3 at < i/ < 1/2.) For instance, at i/ = 1/3, Eqs. (5.7) take the 
form: 

a) C7+^C3<^<C7 + ^C3, 

(5.8) 

16 1 8 

and 

t ^ (^)'' '^-^^ 

in both cases. Provided Eqs. (5.8) are fulfilled, rj respects 1/2 < 7; < 1, and the parameter 
Zc (5.9) can formally acquire any real positive value. 

Therefore, the requirement r^ — Q turns out to be highly interesting since it gives us 
a formal expression for the radius of the domain of localization of the gauge dislocation's 
density profile (just under our approximation (3.29)). 

Now let us focus at the coefficients ri, rs, rg. Here it is appropriate to impose the 
following constraints (notice that J4 = 14 since P2 = 0): 

ri + d2 = C + Xj + ^2 = , 

rs = Je - n + log(A/'')(J5 - %) - = 0, 

re + - J7 - F7 + logA/'(J6 - %) + log'AA(J5 - %) 

+ C:^(l-log(|AA))+di = 0. 

Equations (5.10) simply express the fact that the terms proportional to logp, p^logp, 
and are absent in / (5.1) at p <^ 1 (see expansion (4.19)). Therefore, under our 



(5.10) 



24 



conventions all the terms in (5.2.2) vanish. However, possible contributions oc p^log^p 
should be expected provided / is expanded further. 



B) A/'p > 1 

In this limit ^, the 'integral' contribution /jf in / (5.1) is exponentially small. There- 
fore, the asymptotics of / at some pe, Mpe » 1, is as follows: 



/ ~ /bg(pe) + C 



( ^ 



\2Mpe 



1/2 



sin(jVpe - j) . 



(5.11) 



(2) 

Using (5.11), we obtain the boundary condition a pp \p=p^ — 0, i.e., a free surface boundary 
condition, in the form: 



(2) 

a 



pp 



2k 



P=Pe 



Pi logPe 
8^2 



Pi 



Pe Pe 



where 



Y^{^^p) ^ 



dp 



n 



-M-Y^{Np). 



Besides, under the condition Hp » 1 we obtain: 



(2) 

a 



pp 



2k 



Pi J_ 



- 4di + CN^Yq{Np), 



(5.12.1) 



(5.12.2) 



(2) 

where the representation for a pp is seen from (5.12.1) (with p instead of pe). When 

(2) 

(5.12.1) is fulfilled, i.e., at p = pe, the boundary value of a is given by R.H.S. of 

(2) 

(5.12.2) provided p is replaced by pe- It is seen that the boundary value of cr^,^ tends to 
— 4di at Pe — oo- 

Let us note furthermore that the second equation in (5.10), i.e., the condition r5 = 0, 
can also be fulfilled separately since the difference Jg — Yq depends on the coefficient ^4 
(4.16), which, in turn, contains another free parameter (denoted above as const) to be 
adjusted. Provided r5 = is fulfilled by a choice of const, the following three equations, 
the first and the third Eqs. (5.10), and Eq. (5.12.1), can be written together as a single 
3x3 matrix equation as follows: 



/ 


1 


1 


\ 




( ^1 


1 





Ml 




{'ih 


k 


V2p^ 


1 


M2 


) 







(5.13) 



^AT is large but p is not necessarily large; validity of the replacement (3.29) suggests that «; A/" is 
large 
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where 



^1 = —"Iji 

k = - (Jr - Y,) - logJ\f{J, - %) - log'MiJ^ - n) , 



M..^{l- l„g(lA.)) . 



The determinant (= V) of the matrix in L.H.S. of (5.13), 

- (l-log(|^))^ + f^Pey.(A^Pe) + l, 



(5.14) 



is non-zero. Therefore, Eqs. (5.13) can be solved for di, C, and by the Cramer's rule 
we obtain: 

d, = ((1 - M2)l2 + M,ih - h)) , 

d2 = {{2plMi - M2)h - 2pl I2 + ^3) , (5-15) 

C = V-\h + 2pll2 - h). 

It should be noticed that Eqs. (5.14), (5.15) can be simplified at Afpe 3> 1. 

Therefore, all free constants at our disposal, i.e., di, 0^2 (see (5.1)), C (see (4.18)), 
const (see (4.13.2)), Zc — upc (see (3.29)), C and C" (see (3.30)) are fixed. So far, only 
Q' (3.23) remains to be chosen. It is remarquable that the approach developed allows 

to determine the core radius pc — — sts a, function of two elastic moduli of second order 

(say, p, and u) and of two third order elastic moduli C3 and C7 (under our assumption 
about the simplified density profile). For the given choice of the parameters, all the 
contributions in / (5.1) at p <^ 1 up to those oc p^log'^p (apart from the irrelevant 

(2) (2) 

constant tq) are zero. This means that leading non- vanishing contribution to o"pp, a 
should be expected as oc p^ log'^ p at p — ^ 0. The corresponding term can be deduced 

explicitly after a straightforward calculation: the contribution oc p log p in / must be 

(2) (2) 

accounted for in that case. Thus, a pp, a tend to zero in our model at p ^ 0. Besides, 

(2) 

a pp is zero at the external boundary p — p^. 



5.2 The component a 



zz 



Eventually, let us consider the stress component (Tzz — P- For the background contribution, 
we obtain from (3.24): 

(^bg)., = AA, - Y^(^L^ + Q') ■ (5.16) 
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In the classical approach we consider (cTbg)^^ (5.16) si pc < p < Pe, and we determine Q' 
from the requirement 

azzpdpdcp = 0, (5-17) 



Pc 

which is to express that a "mean" value of azz (i-e., (Tzz averaged over cross-section of the 
bulk) is zero. So, in order to determine Q', we take /bg in the form given by (5.1.2), and we 

(1) 

also take into account Eqs. (3.21), (3.24) where 0bg= (~^/27r) logp is used for evaluation 
of ^^f. Then, using (5.16) and (5.17) we find an expression for Q'. Substituting this Q' 
into (5.16) we obtain: 



'SJzz 



2nJ 1-a 



4^2 1 



1 2 , Pe 

~ 2 2 — 

P Pe - Pc pc 



(5.18) 



where i' = a(l — a) and (1 — a) 



l + u. 



Pi 



Using Eqs. (3.8) and (3.9), one can re-express the coefficient k —— (where pi is given 
by (4.10) with c in it given by (3.42), and is introduced in (4.2)) as follows: 



Pi 
4^2 



foy p 



2tx) 2(1 -u) 



n 
Ap 



- Ap\l + u)C^ 



(5.19.1) 



hV( l-2u 2m' + n'^ 
p + 



27r 



1 - V 



(5.19.2) 



Furthemore, let us express C3 by means of (5.19.1) and then substitute it into (5.18) 
Now {cT\^g)zz acquires the form: 



•SJzz 



, 4^2 



h \2n'' 



1 _ 2 Pe 

2-kJ 8 J \p^ pI- Pc pc 



(5.20.1) 



Pi 

Eventually, we use (5.19.2) to re-express k -^-^ in (5.20), and the stress component (abg 
appears in the form suggested by 0: 



SJ zz 



27r. 



.^ + 4(r^H^-2^)-^9] 



2 , p. 



pt-pc 



PcJ 

(5.20.2) 



It is also important to notice the fact that the integral 



J Nhgpdp 



(5.21) 



pc 



is zero, since the choice of di, d2 in /bg (5.2) ensures vanishing of dfi,g/dp at p = pe,Pc 

(2) 

(the boundary condition for a pp at p = Pe,Pc)- Therefore, (crbg)^^ ^^^^ t)e written in 
the following form: 



Kg), 



-z/A/bg - 2p(l + z/) 



zz 



pI- pI 



"^llpdp 



(5.22) 
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Equation (5.22) is completely equivalent to (5.18). In view of (5.22), it is more clear, how 
Eq. (5.17) is fulfilled. 

Now let us turn to a^z of the modified defect given by (3.35): 

^ -uAf - 2ii (1 + i/)(*33 + ^) + «2 (1 + u) if - /bg) . (5.23) 

(2) (2) 

Let us consider the asymptotical properties of azz- Using the expressions for a pp, a 

(1) 

(5.2.2), the expansion (4.19), and (f) to express ^33 (as well as the nullification conditions 
(5.10)), we obtain that behaves as oc (ri logp+r2 log^ p) plus the contribution tending 
to zero at p — > 0. At M p 3> 1 we obtain: 



1-1/ 



V 



. (5.24) 



(2) 



With the help of (5.12.2) and (5.24) we obtain the trace of «t (and thus the total trace, 
since cr for the screw dislocation is traceless) at Afp » 1: 



tr(<T) 



hCa) = -4(1 + u)di 



+ 



1 + 
1-u 



2t:J 



2p? 



(5.25) 



Besides, C and di given by (5.15) must be substituted into (5.24) and (5.25). Equations 



(2) 



(3.6) and (3.21) demonstrate us how an analogous estimation for the trace of e can be 
deduced. Indeed, 



T 

A e 



(3Ci + C^)h(S) + h{;£) 



= (3Ci + C4) US) - p? (3C3 + C7) {dp 0)' . 



(5.26) 



(1) (2) 

However, at K,p ^ 1, wc estimate {dp (p)'^ b"^/ (27rp)^, and the trace of e can be deduced 
by means of (5.25) and (5.26). 
Eventually, using 



+ 



- 1 h0,Pc]iP), 



(2) 



and the fact that pAf integrated over p from to Pe vanishes, we determine Q'. Sub- 
stituting Q' into (5.23), we find: 

^ 2 

*33 2 j ^33Pdp 



i/A/ - 2p{l + v) 



- 2p{l + u) 



+ 



Wn^pl Sn^pl 
2 

Pe 



^ - 1) ^[o,Pc](p) (5.27) 



2 f 

/ - /bg 2 ~ fhg)pdp 
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We shall not elaborate this expression further. It is enough to notice that the integrals in 
(5.27) are convergent at the lower bands, and thus a^z (5.27) averaged over cross-section 
of the bulk is zero. 

Before to conclude, let us briefly note another possibility which concerns the choice of 
the parameters. Namely, the requirement rg + (ii = can be left aside. In this case, the 

(2) (2) 

stresses cXpp, a (5.2.2) will demonstrate a tending, at p — 0, to the constant values 
~ (re + di) (see (5.2.2)). However, in this case it can be assumed that C = 0. Then, 
instead of (5.13), we shall get just two equations to determine di and d2. First, we obtain 
d2 = h = —ij = —Ti. Then, Eq. (5.12.1) takes the form: 



and it gives us 



h do 

4 - ^ - 2di = , 

Pe Pe 



, ^ ^3 - d2 ^ Pi \0gpe 

'~ 2pl - pI + 2pl 



(2) (2) 

In this case, the asymptotical behaviour of cr^p, cx^,^ is missing the unconventional con- 
tribution due to Yq. Let us stress again that the parameters di and d2 are still different 
in comparison with the analogous conventional results. 



6 Discussion 

A model of non-singular screw dislocation lying along an infinitely long cylindric body 
is investigated in the present paper in the framework of three-dimensional r(3)-gauge 
approach |^. The gauge part of the total Lagrangian is chosen in the Hilbert-Einstein 
form, while the elastic contribution to it corresponds to the energy of elastically isotropic 
continuum given by the terms of second and third orders in the strain components. In 
other words, a second order elasticity approach is adopted in the present paper. 

As it was noticed in P, second order consideration in the framework of the model PP 
would merit attention as an attempt to clarify perspectives of such rather non-traditional 
approach to defects in solids as the gauge Lagrangian approach (based, for instance, on the 
groups either T(3) or IS0{3) = T(3)E) 5*0(3)). Elaboration of related technical details 
could clarify the gauge strategy itself concerning a choice of Lagrangian's constituents, 
of dimensionality of the specific problems, of resolving ansatz, etc. On the other hand, 
it is also intriguing to use such a widely acknowledged and fruitful method as the stress 
function approach [2^ , j3j within an unconventional non- linear gauge framework. Although 
the available gauge solutions of the first order CO, [IDl, [U], [21], [HI] seem to be 
promising, second order consideration could open new aspects of the problem of the 
gauge description of dislocations. 

Let us remind that the linear approach developed in [11 leads to the modifed defects 
which are characterized by the fact that singularities of the ordinary straight dislocations 
are smoothed out. After the classical attempts [2], [3], [1], [S], higher corrections to 
the law 1/p of linear elasticity are known. However, a cut-off near the dislocation axis 
inevitably occures [S], [1]. As it is shown above, the gauge approach P allows to extend 
the description of static screw dislocation to the whole cylindric body containing the 
defect. A use of an approximated density-profile comes to play, and an expression for the 
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radius of the domain of localization of the defect's density profile by means of the second 
and third order elastic moduli appears in the picture proposed in the present paper. 

Second order consideration developed above allows to avoid a stress-free boundary 
condition at an inner radius corresponding to the core radius p = pc- Besides, it allows to 
fix the radius p = pc as a function of second and third order elastic moduli. It removes a 
cut-off which occures in [3^, and the stress components app, a^^ turn out to be continuable 

(2) 

towards the tube's axis. As in the classical approach, it is necessry to subject cr to a 
free-surface boundary condition at the outer radius p = pe- Thus we obtain the solution 
describing a finite cylinder with nonsingular screw dislocation along its axis. Sufficiently 
far from the core, the analytical form of the gauge stress potential of second order found 
above is rather close to the conventional one 

Two possibilities are pointed out for the choice of the parameters in the solution found. 

(2) (2) 

These possibilities enable a pp, a to tend either to constant values or to zero. In the 

(2) (2) 

first case, the analytical form of a pp, a in the region Pc < p < Pe is the same as in 
jH], but the coefficients are nevertheless different. In the second case, an unconventional 

(2) 

contribution is present. However, a is logarithmically divergent (the classical divergency 
is ~ p~^), at p — > 0, i.e., in the region where the approximated form of the density profile 
is most inadequate. In the last case, a weak three-dimensionality may be of help. Some 
estimations which involve the crystallographic parameters are desirable to make contact 
with the known interpretations of the characteristic length k,~^ in terms of interatomic 
spacing (translational gauging), j2H] (non-local elasticity). 

A gauge approach close to ours has been proposed in the series of papers jTU], [TT] , 
|24j . [3T], which is based on the r(3)-gauge Lagrangian written as a combination of the 
terms quadratic in the torsion components. As to the elastic Lagrangian, it is written 
in ^U], pi], |21j without third order terms (since only the linear problems are studied). 
However, it is proposed in |31j to use the terms in the Lagrangian which are related to 
the energy potential of the rotation gradients. Incorporation of such terms in [3T] allows 
to improve the solution found in pP for that modified defect which demonstrates how the 
singularity inherent to the classical edge dislocation is smoothed out. In the far field, 
the stress components found in [21] correctly reproduce those of the edge dislocation. 
The Hilbert-Einstein Lagrangian is highly suggestive representative among the gauge 
Lagrangians of the differential-geometric origin. It leads to the self-contained pictures for 
the modified defects which avoid the singularities of the convetntional solutions. However, 
the contributions of mechanical origin also merit consideration, and further efforts in this 
direction are also needed. 
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Appendix A 

Appendices A and B provide some intermediate results which are helpful in obtaining the 
final asymptotical expressions for the modified stress potential. 

First of all, we directly obtain expansions for KQ{t)lZ{t) and t^^ Io{t)7l(t) at t ^ 1 



'P2 



+ (p4log't + + Ps) logt + + ... 



t-'Ko{t)n{t) ~ -Pi \ogQt) - (palog't + h logt + fca) t"^ 



■2 



{pilog^t + kslogH + h logt + h^t + 



(Al) 



(A2) 



where the coefficients pi, p2, Ps, Pi, and p^ are given by (4.10), and the coefficients ki, k2, 
k^, ki, ^5, are expressed by means of pi, . . . ,p5 as follows: 



ki = ^ + log(|) P2 + 



P3 



^2 = -- + log(-) - + 



4 



P3 



k3 = ^ + log(|)p4 + P5, (A3) 



4 

ki = k - (l - log|) ^ + log(|)p5 



Now let us obtain estimations for the integrals in (4.13). Using (Al) and (A2) we 
obtain at 2; <^ 1: 



f dt 

B{s) - j /o(^)7^(^) j s=.o lo + ^iz^^ + I2 log' 2; + T3 log^ 

s 

+ (J4 log' z + J5 log 2; + Je) 2;' + . . 



where 



Jo = const - Y , ^ ^ Y ' ' ^ ' ^ ^ 

2 ' ^ 8 



(A4) 



X2 = 


-P2 


, X3 = "^^ 


2 


4 


5 


2X6 


= -k - I5; 



(A5) 



and 



Ko{t)n{t) y ~ /Co + + /C22-' logz 

z 

+ /C3 log'z + /C4 log' 2 + /C5 logz (A6) 

+ (/C6 log^ 2; + /C7 log' Z + /Cs log 2; + /Cg) , 
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where 



/Co = T;^ + (l + 2 1og^)^, /Ci = -(l + 21og 



2' 4 



7\ Pi 
V 4 ' 



3 /C3 -fcs + ^4 3 
= -^P4 + y , /Cs = ^ + ^P4, 

h - h 3 
/C9 = ^^ + y-gP4. 

Besides, the coefficient X^- in /Cq is given by the regularized value of the integral: 

00 dt ^ I 

Ik = I Ko{t)n{t) J + I Koit)n{t) +p,\ogQt)t-^ 



(A7) 







dt 

T 



(A8) 



+ k2 + h log t + p2 log^ t 
Using expansions (A4) and (A6), we obtain the expansions we are interested in: 



z 

B(o) - I io(t)n(t) 



dt 

T 



+ % log^ z + K4 log^ z + % logz 
+ (^6 log^ z + % log^ z + Ks\ogz + % 



9 Z 



(A9) 



where 














-log 


(i)Xo + 




i^i = - 


log(|)Xi, K2 = -I, 






K4 = 


-l0g(|)X2 - 




= -Xo-|-log(|)X3, 


% = 


4 


— Z4 , 


= (1 - lo^ 


'2>' 4 4 


- log(|) 2^4 - Xg , 




4 


64 ^ 




-logg)X5- 




% = 


(1- 




+ (|-logi) 


64 ^V2^ 





Analogously, 

00 

/dt ^ 



(AlO) 



2; ^ logz 



+ h log^ -2 + /4 log^ Z + 15 log Z 
+ (/g log^' ;2 + /7 log=^ ;2 + /g log 2; + Tg) 
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(All) 



(A12) 



where ^ 
To = }Co + ^ , 

^^^T + 64 ^^'^ 
Eventually, we obtain the coefficients characterizing the asymptotical behaviour of the 
combination 



h 


— ^1 




4 


+ ^5 


4 = ^ + K 


h 


_ /C2 
64 





(^) y XoW^W y + i^o(^) S(0) - j Io{t)n{t) 



dt 
T 



by summing up the corresponding coefficients given by (A9) and (All). We obtain that 
K2- T2 = (see (AlO) and (A5) for K2, as well as (A12) and (A7) for /a), and the 
coefficients go, ■ ■ ■ , ?7 (4-15), (4.16) appear as follows: 

Qi ^ Ki - Ti , at i = 0, 1 ; 

Qi = Ki+i - Ti+i, at i = 2, ...8. 

Thus, the solution G{z) given by (4.13), (4.14) is estimated, and the final answer is given 
by (4.15), (4.16). 

Appendix B 

First of all, we obtain expansions for tJo{N't)G{Kt) and tYo{N't)G{Kt) at t -C 1: 



k-^ t JoiAft)G{Kt) + q^tlog^iAft) 



+ (ni \og^{Mt) + 712 log(7Vi) + ria) t 

+ (n4 \og\Mt) + 715 log'(AAt) + (...) \og{Mt) + (. . .)) e + 



(Bl) 



.7 



+ (^2 log^(A^i) + mi log^(A^t) + m2 \og^{Ut) + log(A^t) + t 
+ (ms log^(A/'i) + me log^(A/'i) + . . .) + . . . , 



(B2) 
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where ^ 

ni = 3 log(--=) 92 + ?3, 



712 = 3 log^(;^) ?2 + 2 log(-^)g3 + 54, 

A/"^ gi iS/'^X ,2/'*\ ,/^\ 

""^^ ^ ^° " ^ T + W>> + + ' (B3) 

»/-2 92 , 2 

714 = - yV ^ + 95 , 

715 = - jV^ J + 96 + 3 iog(^) «4 , 



and 



mi = log(^) 92 + rii , 7712 = log(^) JT-i + "-2 , 
"^3 = log(|) "-2 + r^s , '^^ ^ ^ ^^^(i) ' ^^^^ 

77l5 = 7l4 , me = A/"^ ^ + log(^) 7l4 + 7l5 . 

Using (B1)-(B4), we pass to the estimation of the integrals which enter into Ij^{p) 
(4.18). First, we obtain: 

oo 

J Jo{J^t)G{Kt)tdt ~ Jo + Ji log(AAp) 
p 

+ {J2 \og\Mp) + Ja \og\Mp) + Ja \og{Mp) + Js) 

+ (Je log=^(AAp) + Jr \og\Up) + . . .) p^ 



where 



2 

^ ,/392 nl^ n- 1 ( ^92 , rii - 7i2 ^ 

^392 7l2 - ni 7l3\ „ ,714 



(B6) 



i4 7*5 ^ 

16" ^ T 

Besides, dots in (B5) imply terms proportional to the first and zeroth powers of \og{Np). 
The constant Xj which gives Jq will be presented below. Further, we obtain: 

00 

j %{m) G{Kt)tdt ^yo + yi WiUp) + y2 \og{Up) 

+ {y^ \og\Np) + y^ \og\Np) + y,, \og\Np) + y^ \og{Np) + y^) p" (^'^) 
+ (ys \og\Mp) + 3^9 \og\Mp) + . . .) p^ 
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where 



3^a = Kf -^ + ^^), (B8) 

and dots in (B7) corresponds to the second, first, and zeroth powers of \og{Hp). 

With the expansions (B5) and (B7) at hands, we pass to the products we are interested 
in to express Ij: (4.18): 

oo 

Yo{Np) j Jo{Nt) G{Kt)tdt ^Yo + Yi log\Afp) + % \og{N p) 
p 

+ (fa \og\Hp) + % \oi{Hp) + % \oi{Hp) + % \og{Hp) + %) p^ 

+ (fs \og\Mp) + % \og\Mp) + . . .) p^ 

where 



(B9) 



% = log(|) Jj, Y^ = Ji, ^2 = Jo + log(|) Ji, % = J2, 

%^J3+ logg) ^2, n = Jl + l0g(|) Js + J4, 

n = ^ (- Jo + (l - log|) Ji) + log(|) J4 + J5: (BIO) 
Yr^—{l- log|)jo + log(^) J5, = -^^2 + Je, 

^9 = ^(1 - log^)j2 - ^ J3 + logg) Je + Jr. 

Now we obtain the following expansions: 

00 

JoiAfp) J YoiAft) G{Kt)tdt = Jo + Ji log'(A/'p) + J2 log(A/'p) 
p 

+ (J3 log^(A/'p) + J4 log3(A/'p) + J5 log2(A/'p) + Je log(A/'p) + J7) 
+ (Jg log^(AAp) + Jg log3(AAp) + . . .) p^ 

(Bll) 
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where 

Jo = Iy, Ji^yi, 1,2,3,4:, 



(B12) 



where 



J7 = -^yo + y?. 

Now we sum up expansions (B9) and (Bll): 

00 00 
/^(p) = - Yo{J\fp) J Jo{Mt)G{Kt)tdt + Jo{Np) J Yo{Nt)G{Kt)tdt = 

p p 
= Jo - fo + {Ji - Y^) log\Np) + (J2 - ^2) log{Np) 

+ (( J4 - %) \og\Np) + ( J5 - fs) log'(AAp) + ( Je - fe) log(AAp) + J7 - %) p" 
+ ((J9 - fg) \og\Np) + 

(B13) 

Jo-yo = :z:y-iog(|)Tj, 

Je - Fe = ^ (jo - (1 - log|) Ji - 3^2) - log(|) J,- 
-J5 + ye, 

Jr-Yy^ (y^ + (1 - log|)jo) - log(|) J5 + y7. 

The terms corresponding to J3 — 13 and Jg — ?8 do not appear since 3^3 = J2 and = n^. 
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Eventually, it is necessary to re-arrange the series (B13) as follows: 
Ir{p) ^Xy- \og{^N)Xj + k^Xog^p - Xj logp + k^pHoip^ 
+ ((J4 - n) \og{N^) + J5 - %)pnoip 

+ ((J4 - F4) 3 log2 + ( J5 - F5) log(Ar2) + Je - %) p2 logp (B15) 
+ ((J4 - log'A^ + (J5 - log'A^ + (Je - ^e) logA^ + J7 - ^7) 
+ (J9 - fgjp'logV, 

where Jj — 1^, i = 4, . . . , 9 are given by (B14) (where (B6) and (B8) must be used) and 
the following notations are adopted: 



Xj = Xj + k%\ogN , 



and 



(00 1 
J Jo{Xt) G{Kt) tdt - jj 
e e 

foo 1 
/ %{m)G{Kt)tdt - k% ( \ogQ-J\ft)- 



(B16) 



(B17) 
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